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Abstract. The theory of singular dislocations is placed within the 
framework of the theory of continuous dislocations using de Rham cur- 
^ , rents. For a general n-dimensional manifold, an (n— l)-current describes 

■ a local layering structure and its boundary in the sense of currents rep- 
resents the structure of the dislocations. Frank's rules for dislocations 

■ follow naturally from the nilpotency of the boundary operator. 

i— i 



1. Introduction 



The aim of this work is to establish a precise relationship between the 
theory of continuous distributions of defects and its discrete counterpart. 
Historically, the latter was developed first by pioneers like Volterra and 
Somigliana. The continuous theory was arrived at later by, among others, 
Bilby |BBS55| . Kroner |Kro59j . Kondo |Kon55| and Noll |Nol67j . Method- 
ologically speaking, the passage from the discrete to the continuous theory 
was perhaps spurred by the realization that certain differential geometric 
objects already provide a heuristic path to generalize the discrete, intu- 
itively graspable, picture. The clearest example is provided by the lack of 
| closure of a Burgers' circuit enclosing an edge dislocation in two dimensions. 

The picture of this event so much resembles that of the lack of commuta- 
tivity of two vector fields, that one would be remiss to ignore the analogy. 
And, in fact, the analogy is in this case fully justifiable on physical grounds. 
In the infinitesimal limit, the lack of closure alluded to above is a lack of 
integrability, whose intensity is measured by the torsion of a distant par- 
allelism associated with the smeared-out underlying crystal lattice. Once 
this particular instance was exploited, the temptation could not be resisted 
to attribute some putative physical meaning to all kinds of other geometric 
objects, from Riemann curvatures to Einstein tensors. The works of Kondo 
and, later, Noll inaugurated the emergence of the opposite paradigm. In- 
stead of building the theory, as it were, from the bottom up, they adopted 
the puristic tenet that the presence of defects (or inhomogeneities) in a con- 
tinuum should be encoded automatically within the constitutive equations 



2000 Mathematics Subject Classification. 74A05; 74E20; 58A25. 

Key words and phrases. Continuum mechanics; dislocations; differential forms; singu- 
larities; de Rham currents; Frank's rules. 

1 



GEOMETRIC ASPECTS OF SINGULAR DISLOCATIONS 



2 



of the body and only there, without any spurious intervention of an atomic 
substrate. The two approaches do not necessarily lead to the same results. 

Since both approaches mentioned above have undeniable merits, it would 
be futile to argue for one to the exclusion of the other. Our intention is, 
therefore, to place the theory of singular dislocations, i.e., dislocations that 
are concentrated rather than distributed continuously, within the elegant 
framework of the smooth theory and so to come full circle to the original 
physical picture. The mathematical tool for this purpose consists of a weak 
reinterpretation of the differential geometric objects of the smooth theory 
in terms of the geometric theory of de Rham currents [dR84, Fcd69j. 

The framework presented below applies to general manifolds. In partic- 
ular, no Riemannian structure is assumed. The theory of integration of 
differential forms on manifolds is used throughout and the necessary back- 
ground may be found in EpslO Segl2| within the context of continuum 
mechanics and in standard references on differential geometry, such as those 
cited therein. Specifically, for the smooth case, we consider a local layer- 
ing of the material as represented geometrically by a 1-form <p, the layering 
form, on a manifold <M which may represent a material body or its image 
under a configuration in space. The condition that the body contains no 
dislocations is thus represented mathematically by the local integrability 
condition dip = 0. This condition implies that at least locally (and also 
globally if ^# is contractible to a point) <p is the gradient of a function u. 
The layers, i.e., hypersurfaces of constant values of u, may be thought of as 
deformed Bravais planes of a crystalline body. Dislocations are present in 
regions where d(p 7^ so that such a system of layers is not available. 

In the non-smooth situation, we generalize the layering 1-form to a de 
Rham (n — l)-current T and the condition for the body to have no dis- 
locations is generalized to dT = 0, i.e., that the boundary of the current 
vanishes. It is noted that in his exposition on singularities in the deforma- 
tions of solids, Cermelli |Cer99] makes use of de Rham currents. 

It is interesting to note that the language of differential forms and cur- 
rents rather than that of frame fields, enables one to analyze dislocations 
associated with a single layering form or current (rather than three). From 
the physical point of view this means that a single family of Bravais planes 
is sufficient for the study of the possible presence of dislocations. It is also 
noteworthy that Frank's conservation rulefl follows naturally from a basic 
mathematical property of currents, specifically, the vanishing of the bound- 
ary of a boundary. 



*In the words of |Rea53[ p. 34], Frank's rules are: 1. "The Burgers vector is conserved 
along a dislocation"; 2. "The sum of the Burgers vectors of the dislocations meeting at a 
node must vanish." 
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2. The smooth theory of crystal dislocations 

2.1. Frames and coframes. The simplest theory of continuous distribu- 
tions of dislocations assumes that a material body 38 is endowed with a 
distinguished smooth field of bases of its tangent spaces. In general, for an 
n-dimensional manifold such a field can be regarded as a section a of 
the frame bundle F^tf, namely 

a : Jt — > FJt 

x i — > a{x) = {ei,...,e„}, (2.1) 

with 7T o a = id^, where tt is the bundle projection and id^ is the identity 
map in We remark that smooth, or even continuous, global sections may 
not exist in general. A manifold for which a smooth section of F^ does 
exist is said to be parallelizable, and the field of basis induces a teleparallelism 
(or distant parallelism) in the manifold. 

There is an alternative (dual) way to look at a distant parallelism. In- 
deed, any frame (that is, any basis of the tangent space T x ^) induces a 
unique dual basis (or coframe) of the cotangent space and vice versa. 

Therefore, a distant parallelism can also be seen as the choice of a particu- 
lar coframe field. Put differently, a distant parallelism induces an M n -valued 
one-form on ^ . Moreover, this coframe field can be regarded as a cross 
section a* of the coframe bundle of 

Let x l (i = 1, n) be a coordinate system. Then, the frame field is given 
by: 

■ d 

e « = e «&^' a = l,...,n, (2.2) 

where the matrix with entries e i a = e i a (x 1 , ...x n ) is nowhere singular. Ac- 
cordingly, the corresponding coframe is given by 

e a = e a i dx i , a = l,...,n, (2.3) 

where e" are the entries of the inverse matrix. As already pointed out, a 
coframe is clearly represented as a collection of n pointwise linearly inde- 
pendent 1-forms indexed by a. 

2.2. The interpretation of covectors as Bravais hyperplanes. Cov- 
ectors may be used to describe geometrically collections of hyperplanes. Let 
W be a finite dimensional vector space and let / € W* be a covector, that 
is, a linear functional / : W — > R. Then, / may be identified uniquely with 
the collection Hji of vectors w in W such that f(w) = 1. This collection 
of vectors constitutes a hyperplane. The difference between any two ele- 
ments of Hf \ belongs to the kernel Hfo of the operator /. Clearly, parallel 
hyperplanes Hjk will be obtained if we consider the elements w € W such 
that f(w) = k for any integer k. Conversely, for any w € W, f(w) may 
be interpreted as the amount of hyperplanes that the arrow representing w 
penetrates. If a covector is multiplied by a positive number a, the density 
of the planes is multiplied by a. A covector includes a choice of orientation 
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for the various planes (a positive side versus a negative side of a plane) and 
the multiplication of / by — 1 reverses the orientation. 

It is natural therefore, to use a covector as a continuous model of a system 
of parallel planes in a crystal (Bravais) lattice. In fact, if W = M 3 , for a 
covector / given as 

/ = fjdx\ (2.4) 

relative to the standard dual basis dx*, the components (f\, f'2, f'3) are pro- 
portional to the Miller indices for the system of parallel planes. Whereas 
the Miller indices are normalized to provide a direction only, the covector 
/ contains additional information as to the density of the layers. Let n be 
the dimension of W. Then, a collection of n linearly independent covectors 
{f a }, a = 1, . . . ,n, induces a collection of n families of parallel hyperplanes. 
These covectors, therefore, represent a Bravais lattice. 

So far we considered a single vector space W. On a manifold the 
interpretation just described can be applied in a point-wise manner, that 
is, to each tangent space T x ^#. Thus, at each point x € the covector 
representing the "direction" and density of the layers is given by the value 
4>{x) of a 1-form (f> : ^# — > T*^#. Noting that this setting does not require 
any additional structure, metric or otherwise, it is natural to study the 
geometric properties of such a differential form as a representation of the 
structure of a class of layers in a lattice. We will refer to such a form as a 
local layering form. 

2.3. Integrability. 

2.3.1. Intuitive considerations . Let F(x) = F*(x)e, be a vector field in a 3- 
dimensional Euclidean vector space with Cartesian coordinates (a; 1 , a; 2 , a; 3 ), 
where {ej} is an orthonormal basis. We recall that the condition that the 
vector field be conservative, namely, that there exists a scalar potential 
function u(x) such that 

F ' - <«> 

is 

dFi - (2 6) 

dxi dx l 

for all i,j = 1,2,3. The above well-known scheme may be generalized to 
an arbitrary differentiable manifold ^ using the terminology of differential 
forms. In particular, one says that a differentiable r-form <p is exact if 
there is an (r — l)-form a such that 4> = da, where d denotes the exterior 
differentiation of forms. Thus, a is a potential form for <p. Every exact form 
4> is automatically closed, that is dtp = 0, since the d operator enjoys the 
property d 2 = 0. Conversely, if the manifold ^# is contractible to a point, 
every closed form <f> is exact, that is, derives from a potential. In the general 
case, when the manifold is not necessarily contractible to a point, if 4> is 
closed, for each x € ^ there is a neighborhood where <\> is exact. Thus, the 
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condition that the form be closed is a generalization of the condition (|2.6p 
for the existence of a local potential. 

In particular, a 1-form eft on an n-dimensional manifold ^ has n compo- 
nents, just as a vector field. Let 4> be represented locally as 

<j> = ^dx\ (2.7) 

where now x l is a manifold coordinate patch. Its exterior derivative is the 
2-form dej) represented locally by 

o ( ^.i dxi A dx ' ~ yZ §iA xi A dx i ] , ( 2i 



»>j \ hJ hJ 



y 



where commas indicate partial derivatives. Thus, the condition d(j) = is 
the analog 

<Pi,j ~ <f>j,i = (2.9) 

of Equation fl2~B . 

We conclude that a closed 1-form represents a locally coherent collection 
of layers that may be identified with deformed lattice planes. In general, 
the 2-form 5 = dip is a measure of the nature of dislocation density and we 
will refer to it as the dislocation density form. 

Let Z be a 2-dimensional manifold of ^# with boundary Y = dZ. Then, 
by Stokes' theorem 



1= <p= 5 (2.10) 

JY Jz 

It follows that J Y ^ is independent of the particular submanifold Z. If there 
exists a submanifold Zq on which dip = 0, i.e., there are no dislocations on 
Zq, then, J z 5 = for any other submanifold Z with boundary Y, even if Z 
passes through a region where dislocations exist (i.e., 5 ^ 0). In the general 
case where no such manifold Zq exists, j z 5 is still independent of Z, and I 
above is a measure the total dislocation embraced by Y in analogy with the 
Burgers vector. 

2.3.2. Parallelism and coordinate systems. In the definition of our distant 
parallelism, the differentiability of the cross section has played no role what- 
soever. If, on the other hand, a* is differentiable, we may calculate its 
exterior derivative da* . In components, we obtain 

T a = de a = e\j dx j A dx\ a = l,...,n, (2.11) 
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namely, r = da* is an M n -valued two-form which we call the torsion of the 
parallelism^ The identical vanishing of the torsion form, namely, 

t£ = a,i,j = l,...,n, (2.12) 

is necessary and sufficient for the existence of a local coordinate system 
such that the original frame field a becomes its natural base. From the 
point of view of the theory of dislocations in continuous media, if the frames 
represent crystalline bases, the vanishing of the torsion implies that the body 
can be smoothly brought to a configuration in which all crystal bases within 
a coordinate patch are mutually parallel, so that there are no defects in the 
lattice. The body is then locally homogeneous. Conversely, a non-vanishing 
torsion is an indication, and perhaps a measure, of the dislocation density 
(or inhomogeneity) . For comprehensive treatments of the general theory of 
inhomogeneity see e.g., [Wan67, EEOTj. 

In the Bravais-lattice interpretation of Section f2.2l one may ask whether, 
given a 1-form </>, there is a "potential" function u : ^# — > M. such that 
(j) = du, where d denotes the exterior derivative (which is identical to the 
gradient in this situation). Such a potential function, if it exists, will label 
the various layers, at least locally. Indeed, these layers would be precisely the 
(local) level surfaces of this potential. If each one of a collection of n 1-forms 
{4> a }, whose values at each point x G ^# are linearly independent covectors, 
derives from a local potential function u a , the values u a {xo) represent a 
point xq £ ^# uniquely and the body has acquired locally the crystalline 
structure of a perfect (non-dislocated) Bravais lattice. Recalling that the 
condition for the existence of a local potential u a for a 1-form <p a is the 
equality of the cross-derivatives, i.e., 




(2.13) 



we recover the integrability condition (|2.12|) . 

As just indicated, the vanishing of the torsion forms is necessary for 
the integrability (holonomicity, homogeneity) of the coframe field e a . In 
terms of the original frame field e a , on the other hand, it is well known 
that the existence of an adapted coordinate system is guaranteed by the 
commutativity of each pair of base vector fields, namely, 

L ea e p = [e a , e p ] = 0, (2.14) 

where L u v is the Lie derivative of the vector field v in the direction of 
the vector field u and where [u, v] denotes their Lie bracket. In terms of 
components, this can be written as 

dei . de i 

l^zl-W^ 4 = 0. (2.15) 
g xJ dx^ y 1 



^Notice an inessential difference with the usual definition, whereby the form takes values 
in T x jjt rather than in R 3 . 
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Since 

the Lie bracket can be expressed as 



= ~-^&L = 3 e? + i ei, (2.16) 



[e Q ,e^] = ^4je CT . (2.17) 

In the context of dislocations, we may identify the Lie bracket between 
the two base vector fields e a and as the local Burgers vector be- 
tween the corresponding crystal directions. Expressing the Burgers vector 
B a fl in the local basis, we may distinguish between the edge component of 
the dislocation density and its screw component, the former being the part 
contained in the plane spanned by the base vectors e a and e^. 

2.4. Frank's rule as the vanishing of the boundary of a boundary. 

Starting from the notions of simplices and chains in an affine space and their 
generalization to manifolds, one arrives at a fundamental geometrical and 
topological result of clear intuitive meaning. It states that the boundary 
dU of any well-defined domain of integration U must necessarily have a 
vanishing boundary, namely, 

2 U = ddU = 0. (2.18) 

In the theory of smooth differential forms, on the other hand, the oper- 
ation of exterior differentiation enjoys a formally similar property, namely, 
for any differential form uj on a manifold 

d 2 u = ddw = 0. (2.19) 

The relation and consistency between these two identities is mediated by 
Stokes' theorem, 

duj = j uj, (2.20) 



U dU 

where uj is an arbitrary (p — l)-form and U is an arbitrary p-dimensional 
domain of integration. 

We will presently show that a smooth version of Frank's rule for dislo- 
cation branching, [Fra51] (see also [Rca53]), can be obtained as a direct 
consequence of these purely geometric identities. For dimension 3, we ob- 
serve that, since for any given smooth coframe field e a the torsion r a = de a 
consists of 3 exact 2-forms, the integral of the torsion over the boundary of 
any 3-dimensional domain of integration U vanishes, 

r a = 0. (2.21) 

dU 

In the physical interpretation, this implies that there are no isolated dislo- 
cation sources, not even as a smoothed-out approximation. In particular, 
consider a tubular domain U and an arbitrarily small neighborhood V of U 
such that the torsion vanishes in V \ U, then, if we intercept U transversely 
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by means of two (oppositely oriented) lids Si and £2 giving rise to a finite 
tube U, we obtain 



In other words, the integral of the torsion over any tube cross section is 
constant. The same reasoning can be applied to a tube with branches, thus 
providing the smooth version of Frank's rule. We note that the restriction of 
the coframe field to V \ U consists of three closed 1-forms, by construction. 
But, since any curve surrounding the tube is not contractible to a point, 
these 1-forms are not necessarily exact. In physical terms, the body minus 
the tube is only locally homogeneous. Moreover, the integral along any such 
non-contractible curve of the coframe 1-forms gives rise to three constants, 
each of them exactly equal to the integral of the corresponding a-component 
of the torsion over any cross section. 

From the heuristic point of view, as the diameter of the tube shrinks, 
we may impose the condition that the torsion increases proportionately so 
as to keep its integral over the cross section constant and thus recover the 
classical form of Frank's rule. The rigorous mathematical treatment of this 
limiting process will be handled in the sequel using the language of currents. 



So as to generalize the notions just introduced, we define a weak telepar- 
allelism p as an IR n -valued current in the sense of de Rham [dR84[ IFed69j . 
We recall that a de Rham r-current is a linear functional on the vector space 
of C°° differential r-forms with compact supports in such that T(0) — > 
if the components of <j> and all their derivatives tend to zero uniformly in the 
support of 4>. A de Rham current is the natural generalization of a Schwartz 
distribution to manifolds where the forms <j> are analogous to test functions. 
As such, currents provide a tool for the study of non-smooth, concentrated, 
physical phenomena, e.g., dislocation lines and slip surfaces. We identify 
naturally an n-tuple of r-currents with an M™-valued r-current. 

3.1. The current induced by a form. A 1-form <p on ^# may be paired 
with a smooth (n — l)-form ip, having a compact support, to produce a real 
number in the form 



Here, (p/\ifi denotes the exterior product of the two forms, an n-form having 
a compact support which may be integrated over the n-dimensional manifold 
jffl . Thus, the form ip induces a linear functional T v acting on the vector 
space of (n — l)-forms of compact supports in j$ in the form 





3. The weak counterpart 
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If all the derivatives of the local representatives of the form -0 tend uniformly 
to zero in compact subsets of the domains of charts in then T^ip) tends 
to zero. It follows that T v is indeed an (n — l)-current. 

Remark 3.1. The action (|3.2p may be given a physical interpretation in a 
different context. The 1-form ip may be interpreted as a force field per 
unit value of a certain extensive property. For example, as the electric field 
in the case where the property under consideration is the electric charge. 
Thus, the question whether ip is closed corresponds to the question of the 
existence of a potential function for the force field. The (n — l)-form ip is 
interpreted as the flux field of the property under consideration so that for 
any n-dimensional region R C 



is interpreted as the total flux of the property through the boundary dR. 
Thus, Ttpty) in f)3.2f) may be interpreted as the total power expended by the 
force field while the transport of the property is given by the flux field ip. 

Remark 3.2. A collection of n 1-forms in 1", e.g., {e 1 , ...,e n } may be inter- 
preted as the collection of gradients of the components of a velocity field. 
The corresponding velocity field will be incompatible, or contain dislocation 
rates, if the forms are not closed. Thus, for a collection of n (n — l)-forms 
ip 1 the action 



may be interpreted as the mechanical power performed by the stress matrix 
having the components ip 1 (each ip 1 has n components itself) on the velocity 
gradient. 

3.2. Examples of currents in general. An (n — l)-current in the form 
(13. 2p is very special as it is induced by a smooth 1-form. As such it may 
be identified with the form (p. The fine topology used for the test forms 
enables one to define currents which are a lot less regular. Such currents, 
generalizing the local layering forms, will be referred to as local layering 
currents. We present below a number of examples. 

3.2.1. The current induced by a form. We have seen already that a current 
T<p, induced by a closed 1-form ip according to Equation (|3.2p . represents a 
locally coherent system of layers as in Section 12.3.11 . 

3.2.2. Incoherence at an interface. Let j$ = M? = {(x 1 ,^ 2 )} and consider 
the 1-form 



The form (p is not continuous, yet the current T„ as in Equation (|3.2p is 
well defined. For x 2 < 0, tp represents a collection of vertical layers and for 




(3.3) 




(3.4) 




(3.5) 
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x 2 > 0, ip represents a collection of vertical layers that are a times more 
dense. Thus, the form tp describes incoherence at the interface x 2 = for 
a ^ 1. 

3.2.3. A Dirac current. Let v\,... , u n -i be a collection of vectors in T xo ^ 
for a point xq € Then, 

T(ip) = ip(x )(vi,...,v n - 1 ) (3.6) 

is an (n — l)-current which is a generalization of the Dirac delta distribution. 

3.2.4. The current induced by the exterior derivative. Consider the (n — 2)- 
current T dip induced by a 1-form ip as 



T dv (u) = dpAu, (3.7) 

for any compactly supported (n — 2)-form to. Using the basic property of 
exterior differentiation whereby 

d(aA/3) = daA/? + (-lfaAd/3, (3.8) 

for an r-form a, Equation (|3.7p may be written in the form 



T Aip {u) = / d(pAoo) + / pAdu, 

ip Awl / p Aduj, (3.9) 
p> A da;, 



where in the second line we used Stokes's theorem and in the third line we 
used the fact that oj is compactly supported in M . 

3.2.5. Polyhedral chains as currents. A current T s can be uniquely associ- 
ated with an (n — l)-simplex s in It is defined as 

r s (v) = f V, (3.io) 

J s 

for any compactly supported (n — l)-form ip. Thus, rather than a continuous 
system of layers modeled by a form p as in Section f3. 2. II above. T s represents 
a single "concentrated" layer. For example, a simplex s inside jtft may 
represent a cut inside the body where an additional layer of atoms has been 
added or removed. 

Evidently, we may extend this definition to an arbitrary chain A = 
Y^p CLpSp and define Ta by 



T A W = ^2 a p[ ^ (3-11) 

" Sri 
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3.2.6. The product of a current by a function. Let T be a current and u a 
smooth function. Then, one may define the product current uT by 

uT(ip) = T(utfj). (3.12) 

In particular, for the current T s of the previous example, 




(3.13) 



3.3. Dislocations as boundaries of currents. We recall that the bound- 
ary of a p-current T is the (p — l)-current dT defined by 

ar(w) = T(dw). (3.14) 

In case dT = 0, one says that T is closed. Just as a current is a non-smooth 
generalization of the system of layers represented by a 1-form, dT = is a 
generalization of the condition dp = (see Section 13.3.11 below) implying 
coherence of the system. In fact, a theorem by de Rham (see [dR84l pp. 79- 
80]) asserts that any closed current is homologous to a current T v induced 
by some smooth form p, i.e., there exists a current S such that for each 
compactly supported smooth form ip, 

(T-T lp )^)=T(^)- [ ^AV = 3S(V). (3-15) 

In other words, the homological properties of T may be obtained by the 
analogous properties of a current induced by an approximating smooth form 
<p>. As the dislocation structure of a smooth form ip is obtained by dip and as 
8T V = T^ (Section 13.3. 1 p . it is natural to obtain the dislocation structure 
induced by the (n — l)-current T from its boundary dT. Thus we will refer 
to the (n — 2)-current D = dT as the dislocation current. 

Again, we demonstrate the significance of these notions in the following 
examples. 

3.3.1. The boundary of a current induced by a smooth form. Consider Sec- 
tion 13.2.41 above. It follows from Equation (13. 9j) that for a 1-form (p, 

dT v = IV (3.16) 

We conclude that if p is closed then dT^ = 0. This is just the condition for 
coherence phrased in terms of currents. 

3.3.2. The dislocation current for an incoherent interface. Consider Exam- 
ple 2 above. Using M 2 ~ and 1R 2+ to denote the lower and upper half planes 
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in M 2 , we have for each 0-form lo, 

dT v (oj) = dx 1 A do; + / adx 1 A du, 
Jr 2 - Jr 2 + 

d(dx 1 Aw) + / d 2 x 2 A u 



/ ad(dx 1 Auj) + ad 2 x 2 A w, (3.17) 
dx 1 A a; — / adx 1 A a;, 

« 2 " J 

(a — 1) / codx 1 , 



where it is noted that cM 2 and (M 2+ contain the set L = {(a^O)} but 
with opposite orientations. Let Tl be the 0-current in M? defined by 

T L (u) = / udx 1 . (3.18) 



Then, the preceding calculation shows that 

8T V = (a - l)T L . (3.19) 

Indeed, for the case where a = 1, dT v = and is a closed current which 
represents a coherent collection of layers. In case a^l, the dislocations are 
concentrated on the line L which is the support of dT^, i.e., dT v (oj) = for 
any 0-form (a function) uj whose support is disjoint from L. 

3.3.3. The dislocation current induced by a polyhedral chain. Referring to 
Section I3.2.5| we note that by Stokes's theorem, 

dTju) = fdu, 

(3.20) 

u, 

ds 

so that 

8T S = T 9s , (3.21) 

where ds is viewed as a polyhedral chain. Thus, as one would expect, the 
dislocation line is the boundary of the embedded simplex. Evidently, the 
boundary operator is linear and may be extended in this case to a polyhedral 
chain. 

3.3.4. General incoherent interfaces. Let Y be an n-dimensional submani- 
fold of ^# with boundary Z = dY. Let (p be a closed form and consider the 
layering current 

T(if>) = J cup A i\) + J_ (p Aip, (3.22) 
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where a > € R and Y is the manifold with boundary —Z whose interior is 
jffl\Y (so that the orientation of d(^\Y) is the opposite of the orientation 
of Z). We have 

D(uj) = 8T(uj) = I A do; + J if A dw, 

ad(<p A u) + / ad</2 A w, (3.23) 

d(<y? Au)+ I di^Aw, 

and using the assumption that d<y? = 0, it follows that 

8T(u) = (1 -a) / y)Aw. (3.24) 

Thus, the dislocations are distributed over the boundary of Y while the 
material is coherent inside and outside Y. 



Y 



3.3.5. A dislocation line. Consider the case where = (—1, l) 3 is an open 
cube in R 3 . Let 

s = {(0, x 2 ,x 3 ) £ | x 2 < 0} (3.25) 

equipped with the orientation induced by the form dx 2 A dx 3 and let T s be 
the 2-current defined by 

T s (cp)=j V (3.26) 

for any 2-form ip compactly supported in jtft . It follows that for any com- 
pactly supported 1-form u, 



dT s (u) = 




(3.27) 



where L = {(0,0, x 3 ) € ./#} oriented naturally by the form dx 3 . Clearly, L 
represents the line of dislocation associated with the half plane s. Notice 
how closely the layering current T s matches the addition of a half plane of 
atoms as depicted in standard texts on dislocations. 



3.3.6. The boundary of a product of a function and a chain. Using again 
the setting of Section [3.2.61 let s be an (n — l)-simplex in the n-dimensional 
manifold ^# and let u be a smooth function. Set 



T us (ip) = J uip 



(3.28) 
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for every compactly supported (n — l)-form ip on It follows that for 
every compactly supported (n — 2)-form u on ^# one has 

dT us (uj) = T us {Auj), 
= \ udu, 

J s 

= / d(uu) - / du A u, (3.29) 

J s J s 



uuj — I duAw, 
ds J s 

= {T u ds — T s \- du)(w). 

Here, we have used the notation 

Tl ol{uj) = T(ot A oj) (3.30) 

for an r-current T, a p-form a and a compactly supported smooth (r — p)- 
form w. Thus, in general, one has 

dT us = T uds - T s \- du. (3.31) 

Clearly, one may replace s above by a polyhedral chain or a smooth sub- 
manifold (using triangulation) . 

3.3.7. A node of three dislocation lines. Let ^ = (— 1, l) 3 C M 3 , si = 
{(0,x 2 ,x 3 ) G Jl | x 2 < 0, x 3 > 0}, s 2 = {(sS^.O) G ^# | x 1 > 0, x 2 < 0}, 
S3 = {(x^x^O) G ^# I x 1 ^ 0, x 2 ^ 0} where the quarter planes si, S2, S3 
are oriented by the normals n\ = (1, 0, 0), n 2 = (0, 0, 1) and n 3 = (0, 0, —1), 
respectively. Consider the layering described by the current 

T = T aiSl + T a2S2 + T a3S3 , (3.32) 

so that 

T(V) = J> [ V, (3.33) 

i=l Js i 

for a smooth compactly supported 2-form if;. Thus, 
D(u) = 8T{oj), 
3 



i=l J Si 
i=l 



(3.34) 



'dsi 

= ( ai T Ll + a 2 T L2 + a 3 T L3 + (ai - a 2 - a 3 )T L )(w), 

where the one dimensional simplices L p are define as follows: L\ is the 
segment from the origin to (0,0, 1), L 2 is the segment from (1,0,0) to the 
origin, L3 is the segment from (—1,0,0) to the origin, L is the segment 
form (0, —1, 0) to the origin. It is noted immediately that if the dislocation 
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current is supported only on the "fork" L\ U L2 U L3, then, one has the 
condition 

ai = a 2 + a 3 . (3.35) 

Evidently, this result, a particular case of Frank's second rule, will also hold 
if the cube is deformed under any embedding in M 3 . Furthermore, the choice 
of planes is immaterial. (See also Section [3.4,31 ) 

3.4. The boundary of a boundary and Frank's rules. The theory of 
currents provides a generalization of the intuitive result of combinatorial 
topology that the boundary of the boundary of a chain is zero. This follows 
immediately from 

d 2 T(a) = dT(da) = T(d 2 a) = 0. (3.36) 

The dislocation current D is obtained as the boundary of a current T. 
Hence, 

3D = 8 2 T = (3.37) 

is a condition that the dislocation current must satisfy. In other words, the 
dislocation current must be closed. 

We may use this result in the following situations. 

3.4.1. Frank's first rule. Let L be an (n — 2)-dimensional submanifold with- 
out boundary in the manifold (For example, in a three dimensional 
situation, L could be a curve that does not have ends inside We as- 

sume that L is the support of the dislocation current D. We want to examine 
the possibility that the dislocation current is of the form 

D = T uL (3.38) 

for some real valued function u defined on Thus, there is a local layering 
(n — l)-current S such that 

D = T uL = 8S. (3.39) 
Using Equation (|3.3ip for the submanifold L (instead of s), we have 

= dD, 

= dT uL , (3.40) 
= T udL - T l l. du, 

and, by the assumption that dL = 0, we conclude that 

du = (3.41) 

so that u must be constant on L. This result is clearly analogous to the 
Frank's first rule (in which case L is 1-dimensional) and it is an example of 
the constancy theorem of geometric measure theory. 



GEOMETRIC ASPECTS OF SINGULAR DISLOCATIONS 



16 



3.4.2. Frank's first rule for the boundary of a submanifold. Let Z be an 
(n — l)-dimensional submanifold with boundary dZ of let u be a smooth 
function on ^# and consider the (n — l)-current T uZ given as 

T uZ {i>) = ! mp. (3.42) 
Jz 

Using the analog of (|3.3ip for the submanifold Z, we have 

dT uZ = T udz -T z ^du. (3.43) 

Assume that dT uZ is supported on dZ. By applying dT uZ to forms whose 
supports are disjoint from dZ, it follows that du must vanish on Z. We 
conclude therefore that if u is not constant on Z, the support dT uZ contains 
points outside dZ. In the context of dislocations, if the "intensity" of the 
dislocations along a certain line L is not constant, there should be additional 
continuous dislocations on the dislocation surface Z outside L. 

3.4.3. Frank's second rule. Let ^# be a nonempty bounded open subset of 
M 3 and let A 6 jtft . Consider 3 curves Li, i = 1,2,3 in ^ such that each 
Lj is connected in it originates at A and is the intersection of the image 
of a curve Cj : [0, 1] — > R 3 with ^ such that Cj(l) ^ ./#. (In other words, 
each Ci ends on the topological boundary of jtft where it is noted that 
has no boundary as a manifold and not as a current.) Thus, dLi = {A} as 
a manifold and dT^ = Ta (the Dirac delta) as a current. We examine the 
case where the dislocation current is given by D = Y2% a i^Li- It follows that 
for an arbitrary compactly supported smooth 0-form a, 

= dD{a), 




(3.44) 



= y^ai<x(A). 

i 

We conclude therefore that = Yl a i ■ This condition is evidently analogous 
to Frank's second rule for dislocations. 
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